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ABSTRACT 

The Press-Schechter (PS) formalism for mass functions of the gravitationally collapsed objects is reanalyzed. 
It has been suggested by many authors that the PS mass function agrees well with the mass function given by 
A^-body simulations, while many too simple assumptions are contained in the formalism. In order to understand 
why the PS formalism works well we consider the following three effects on the mass function: the filtering effect 
of the density fluctuation field, the peak ansatz in which objects collapse around density maxima, and the spatial 
correlation of density contrasts within a region of a collapsing halo because objects have non-zero finite volume. 
According to the method given by Yano, Nagashima & Gouda (YNG) who used an integral formula proposed by 
Jedamzik, we resolve the missing link in the PS formalism, taking into account the three effects. While YNG 
showed that the effect of the spatial correlation alters the original PS mass function, in this paper, we show that the 
filtering effect almost cancels out the effect of the spatial correlation and that the original PS is almost recovered 
by combining these two effects, particularly in the case of the top-hat filter. On the other hand, as for the peak 
ansatz, the resultant mass functions are changed dramatically and not in agreement especially at low-mass tails in 
the cases of the Gaussian and the sharp fc-space filters. We show that these properties of the mass function can be 
interpreted in terms of the kernel probability P(Mi IM2) in the integral formula qualitatively. 

Subject headings: cosmology: theory — galaxies: formation — galaxies: luminosity function, mass function — 
large-scale structure of universe 



1. INTRODUCTION 

It is one of the most important work to determine the num- 
ber density of collapsed objects from the statistical properties 
of the density fluctuation field at an early stage of cosmic struc- 
ture evolution. Such number density, or the mass function of 
halos, constitutes a main body of analyses of cosmic structures 
such as galaxies. The pioneering work was done by Press & 
Schechter (1974, hereafter PS). They derived a mass function 
by relating a one-point distribution function of density fluctua- 
tion with a multiplicity function of halos. It is very simple and 
useful, and therefore has been applied to many problems related 
with statistics of cosmic objects by many authors. Moreover, in 
order to take into account the formation histories of galaxies, it 
has been extended to estimate the mass function of progenitors 
of a present-day halo (Bower 1991; Bond et al. 1991, hereafter 
BCEK; Lacey & Cole 1993, hereafter LC; Kauffmann & White 
1993; Somerville & Kolatt 1998). Thus to remove uncertainties 
in the PS formalism and to refine it are very important. 

Recent understanding on the cosmological structure forma- 
tion is based on a cold dark matter (CDM) model which is 
a kind of the hierarchical clustering scenario. In the CDM 
model, the Universe is dominated by a cold dark matter grav- 
itationally, and the self-gravitating structures are formed via a 
gravitational growth of small initial density fluctuations. Since 
higher density regions collapses earlier and since the power at 
smaller scales of the fluctuations is stronger than that at larger 
scales, larger objects which are formed by collapse of large- 
scale fluctuations are formed by clustering of smaller objects 
hierarchically. Formation processes of objects such as galax- 
ies and galaxy clusters are qualitatively well understood in the 
context of such hierarchical clustering scenario. However, once 
we would like quantitative estimate, especially in regard with 
the statistical properties such as luminosity function of galaxies 



and abundance of clusters of galaxies, we face on difficulties 
caused by uncertainties to estimate the mass function of cos- 
mological objects. 

The PS formalism is the most useful even at present. This 
formula assumes a Gaussian distribution of density contrast and 
spherically symmetric collapse of objects (briefly reviewed in 
§2). In the PS formalism, however, there is a problem of the so- 
called "fudge-factor-of-two". Because their consideration was 
only in overdense regions, they simply multiplied the factor of 
two in order to consider all matter in the Universe. 

Another approach explaining the distribution of cosmologi- 
cal structures has been developed, that is, the peak formalism. 
This was firstly investigated by Doroshkevich (1970), and ex- 
tensively studied by Peacock & Heavens (1985) and Bardeen et 
al. (1986, hereafter BBKS). In this formula, considering con- 
figurations of local density maxima through the first and second 
derivatives of the density fluctuation, the number density of the 
density peaks is derived. However, peaks in a field with larger 
smoothing scale often include those with smaller smoothing 
scale. Thus the mass function derived from the peak formal- 
ism is in agreement with numerical experiments only at very 
large mass scales (Appel & Jones 1990). In order to derive a 
correct mass function, we must remove such small scale peaks 
contained in larger peaks (the so-called cloud-in-cloud prob- 
lem). By using "the confluent system formalism", Manrique & 
Salvador-Sole (1995) derived a mass function taking into ac- 
count the removal of nested peaks in the case of a Gaussian 
filter. 

Peacock & Heavens (1990; hereafter PH90) and BCEK at- 
tempted to solve the cloud-in-cloud problem by considering a 
trajectory of a density fluctuation against the varying smoothing 
scale. The method is called the excursion set formalism accord- 
ing to BCEK. Assuming the spherical collapse, the threshold 
for the collapse is determined by only the value of the den- 
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sity fluctuation, independent of mass scale. Thus we adopt as 
a mass scale of the density fluctuation of a viriallized object, 
a maximum mass scale at which the density contrast is firstly 
upcrossing the threshold for decreasing mass scale. When the 
sharp fc-space filter is adopted, they obtained the fudge-factor- 
of-two because of the Brownian motion of the trajectories due 
to uncorrelativity of phases of Fourier modes. However in the 
case of general filters, since the correlation between different 
Fourier modes originates the non-Markov behavior, a Monte 
Carlo method or complicated approximations are required, as 
done by PH90 and BCEK. It should be noted fliat this prob- 
lem had already been considered by Epstein (1983; 1984) in 
the case of a Poisson fluctuation of point masses. 

In contrast to the excursion set formalism, which is based on 
a differential formula such as the diffusion equation, Jedamzik 
(1995) proposed another approach based on an integral equa- 
tion. YNG fixed an inconsistency of estimation of probability 
in the original analysis by Jedamzik (1995) and also obtained 
the fudge-factor-of-two in the case of the sharp fc-space filter. 
While within the case of the sharp ^-space filter, they extended 
the Jedamzik' s integral formula to include the effect of the spa- 
tial correlation which should be considered within a collapsing 
region in order to solve the cloud-in-cloud problem. They found 
that the factor of two is lost by the spatial correlation and pre- 
dicted larger number density of halos compared to the PS mass 
function at large scale at which the variance of the density fluc- 
tuation becomes below unity. They also connected the peak for- 
malism with the Jedamzik formula and found that the low-mass 
tail of mass function is changed. Nagashima & Gouda (1997) 
confirmed that the mass function given by YNG is in agree- 
ment with that given by the Merging Cell model proposed by 
Rodrigues & Thomas (1996), in which the density fluctuation 
field is realized by Monte Carlo method and halos are identified 
according to the hnear theory and the spherical coUapse model. 

On file oflier hand, Monaco (1995; 1997a, b;1998) inves- 
tigated the dynamical effect of the collapse by using nonlin- 
ear dynamics such as the Zel'dovich approximation (Zel'dovich 
1970). As well known, the critical density contrast is 1 .69 inde- 
pendent of mass scale for the spherical collapse in the Einstein- 
de Sitter universe. However, the probability finding objects 
collapsing spherically is very small and most of halos collapse 
non- spherically. Taking into account such direction-dependent 
dynamics, if we interpret a first shell-crossing epoch as a col- 
lapsing epoch, it was found that the resultant mass function is 
similar to the PS mass function with effectively decreased crit- 
ical density below 1.69. 

The purpose of this paper is to clarify why the PS mass func- 
tion works well and well reproduces mass functions obtained 
by A'-body simulations (e.g., Sheth & Tormen 1999; Jenkins et 
al. 2001). Hereafter we call this problem the missing link in 
the PS formalism. In this paper we consider the following three 
effects based on the integral formula developed by Jedamzik 
(1995) and YNG. First the filtering effect is investigated. We 
use the sharp fc-space, the Gaussian, and the top-hat filters. Sec- 
ond the peak effect is considered, based on the peak ansatz in 
which objects collapse around density maxima. Third the spa- 
tial correlation is introduced in the same way as YNG, but for 
the above three filters. Finally all the effects are considered si- 
multaneously. Since it is very complicated that the dynamical 
effect as considered by Monaco is introduced simultaneously, 
we will only discuss this effect in §6. 

The structure of this paper is as foUows. First in §2, we de- 



fine the smoothed density field introducing the filters. In §3, we 
review the PS formalism. In §4, we show the formulation by us- 
ing an integral equation, based on Jedamzik (1995) and YNG. 
In §5, we show the resulting mass functions and consider the 
effects of the filters, peaks, and spatial correlation. In §6 we 
discuss some uncertainties which we do not investigate in de- 
tail. We devote §7 to summary and conclusions. Many detailed 
definitions and equations are summarized in Appendix. 

2. SMOOTHED GAUSSIAN RANDOM FIELD 

First, we define the density contrast as S(x) = where 
pix) is the density at the point x and p is the cosmic mean den- 
sity. The Fourier mode of the density contrast, 5^, is obtained 

(1) 



-I 



5k = / 5(x)e'''"dx 



This density contrast is smoothed out by a window function 
Wmir) with a smoothing mass-scale M, 



6m{x) = j Wm{\x' -x\)6{x')dx' 

= (2^)3 / WikR)6r.e-'^-dk, 



(2) 



(3) 



where W(kR) is the Fourier transform of the window function, 
WMir), /? is a scale related to the mass M, and k= |k| . 

The following three filters are often used as the window func- 
tion by using the smoothing length R: 



1. Top-hat filter 



WikR) = 



ikR) 



■ (sin kR-kR cos kR) . 



2. Gaussian filter 



WMir) : 



1 



WikR) = exp 
3. Sharp ^- space filter 



(27r/?2)3/2 
k^R^ 



exp 



(4) 
(5) 

(6) 
(7) 



WMir) : 



sin krr— krr cos krr 



(8) 
(9) 



In^r^ 
WikR) = eikc-k), 

where kc is the cut-off wave number, kc — and Oix) is the 
Heaviside step function. 

Next we define a variance of the density field with a smooth- 
ing mass-scale M, (Jo(M). The variance is obtained by using the 
power spectrum Pik) = (|5kP) as follows, 

1 



a',(M)^{Sij) = 



(2^)3 



WHkR)P{k)ATTk'-dk. (10) 



The one-point distribution function in the Gaussian random 
field is characterized by the above variance. The normaUza- 
tion of the power spectrum is given by a characteristic mass M* 
at which o"o(A^*) = 1- 

In connection with the definition of kc, here we define the 
mass M with a scalelength R or kc. In the case of the top- 
hat filter, the mass M is clearly related to 4nR^p/3. On the 
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other hand, for the other fihers, the relation is not trivial. For 
example, LC proposed that the mass should be defined as the 
enclosed mass by the window function Wtnir) renormalized to 
WM(r = 0) = 1 . It should be noted that this uncertainty does not 
emerge in the case of the Gaussian filter because it can be renor- 
malized into M* in the case of the scale-free power spectrum. 
As for the sharp fe-space filter, in this paper, we adopt a simple 
relation, kcR = n, while LC adopted kcR = (97r /2) '/^ . The kc-R 
relation will be discussed in §6. 

By using the spherically symmetric collapse model (Tomita 
1969; Gunn & Gott 1972), the critical density contrast 6c is 
equal to 1.69 in the Einstein-de Sitter universe. So the region 
with (5 > 5c is already included in a collapsed region. We adopt 
this value throughout this paper. Note that Monaco (1995) 
found that effective critical density contrast taking into account 
a non-spherical collapse is approximated as (5c — 1 .5 by deaUng 
with 6c as a free parameter and fitting the PS mass function to 
that given by him. 

In the foUowings, we restrict ourselves to the Einstein-de Sit- 
ter cosmological model and scale-free power spectra, Pik) oc 
with n = and n = -2, except for §§5.4 in which we wiU show 
mass functions in a CDM universe. 

3. THE PRESS-SCHECHTER FORMALISM 

In this section, we review the Press-Schechter formalism 
briefly. The probabiUty of finding the region whose linear den- 
sity contrast smoothed on the mass scale M, 6m, is greater than 
or equal to 6c is assumed to be expressed by the Gaussian dis- 
tribution given by 



/(> 6c;M) = 



1 



exp 



2a2(M), 



d6. 



(11) 



/InaoiM) J 5, 

This probability corresponds to the volume fraction of the re- 
gion with more than or equal to 6c in the density field with 

smoothing scale M to the total volume in a fair sample of the 
Universe. Therefore, the difference between /(> 6c', M) and 
/(> 6c;M + dM) represents the volume fraction of the region 
for which 6m = 6c precisely. The density contrast of an isolated 
collapsed object must be precisely equal to 6c because an object 
with 6 > 6c would be eventually counted as an object of larger 
mass scale. The volume of each object with mass scale M is 
M/ p. Then we obtain the following relation. 



Mn(M) 



dM = - 



df(> 6c;M) 



dM, 



(12) 



p dM 
where n(M) means the number density of objects with mass M, 
that is, the mass function. However, the underdense regions 
are not considered in the above equation. Hence, Press and 
Schechter simply multiply the number density by a factor of 
two, 

Mn(M),„ ^ dfi>6c-M) ^^^ 

dM = -2 dM. (13) 

p dM 

This factor of two has long been noted as a weak point in the 
PS formula (the so-called cloud-in-cloud problem). PH90 and 
BCEK proposed a solution to this problem by taking account 
the probability P^p that a subsequent filtering of larger scales 
results in having 6 > 6c at some point, even when at smaller fil- 
ters, 6 < 6c at the same point (the excursion set formalism). In 
the sharp fc-space filter, because of the uncorrelation between 
different Fourier modes 6\i, the trajectory of 6m behaves as a 
random walk for adding a new Fourier mode. In this case, they 
found that the factor of two in the PS formula is reproduced. 

' In YNG, we used a notation P(Mi,M2) according to Jedamzik (1995). In this 
probability for i5i above 5c in a Gaussian field Mi within a collapsing region with 



4. INTEGRAL EQUATION APPROACH 

Jedamzik (1995) proposed another approach to the cloud-in- 
cloud problem based on an integral equation. 

Now, we consider the regions whose smoothed Unear density 
contrasts with smoothing mass scale Mi are above 6c. Each re- 
gion must be included in an isolated collapsed object with mass 
M2 larger than or equal to Mi . Therefore, we obtain the follow- 
ing equation. 



fi>6c;Mi): 



f 



M7 

P{Mx\M2)^n{M2)dM2, 

P 



(14) 



where P{M\ IM2) means the conditional probability of finding a 
region of mass scale Mi in which 6mi is greater than or equal 
to 6c, provided it is included in an isolated overdense region of 
mass scale M2.' Hereafter we call the procedure in which the 
mass functions are estimated by solving eq.(14) the Jedamzik 
formahsm. If f(Mi|M2) is simply given by the conditional 
probabiUty p{5m^ > 6c\6m2 = 6c), P{M\ IM2) is written as follows 
by using the Bayes' theorem, 

P(Mi|M2) = p{6Mi > 6c\6m, = 6c) 

1 r r \ {6Mi-6cf \ 



'Ml 



1 

2' 



(15) 



where = (Jq{Mi) - (Jq{M2) by using the sharp ^-space filter. 
Thus we can obtain the PS formula, naturally including the fac- 
tor of two as can be seen from eqs.(14) and (15). The reason 
originating the factor of two is a behavior of the random walk, 
and is essentiiiUy the same as the excursion set formalism. 

However, it is insufficient to use eq.(15) for more realistic es- 
timation of the mass function because it is necessary to consider 
the spatial correlation of the density fluctuations due to the finite 
size of objects. Besides the sharp A:-space filter is used which 
leads (7^y(, to crQ(Mi) - al(M2), then we should evaluate how the 
mass function is changed by other filters. Therefore, we must 
consider the probability P{r,M\ IM2) of finding > (5c at a dis- 
tance r from the center of an isolated object of mass scale M2 
with general filters. Then we can get the probabiUty P{Mi IM2) 
by spatially averaging P(r,Mi|M2). 

Because we expect that the isolated collapsed objects are 
formed around density peaks, the constraints to obtain the 
above probabiUty P(r,Mi IM2) are given as follows: 

1. The linear density contrast, 6m2, of the larger mass scale 
M2 should be equal to 6c = 1.69 at the center of the 
object. 

2. Each object of the mass scale M2 must contain a 
maximum peak of the density field, i.e., the first 
derivative of the density contrast V^Mj must be equal 
to and each diagonal component of the diagonalized 
Hessian matrix C, of the second derivatives must be less 
than at the center of the object. 

3. The density contrast of the smaller mass scale 
Mi(< M2) which has been already collapsed and which 
is included in an object of mass scale M2 must satisfy 
the condition (5m, > 6c at distance r from the center of 
the larger object. 

paper, we adopt a form P(Mi IM2) in order to stress that this shows a conditional 
M2. 
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In order to evaluate the above conditions, in this paper, we in- 
vestigate the following effects: (1) filters, (2) density peaks, and 
(3) spatial correlation. Some other effects wiU be discussed in 
§6. 

5. RESULTS 

5.1. Filtering effect 

In general filters except for the sharp ^-space filter, the de- 
scription of the Markovian random walk in the excursion set 
formalism cannot be adopted. PH90 and BCEK calculated the 
mass function by Monte Carlo method realizing Fourier modes 
Sk in such a non-Markov case. On the other hand, we can eas- 
ily calculate the mass function by substituting the conditional 
probability in eq. (15) to general cases as a simple approxima- 
tion. It is adequate approximation for the purpose of this paper 
in order to understand the various effects mentioned above, at 
least, qualitatively. 

For this purpose, what we do is to estimate the probability, 

P(Mi IM2) = pi5M, > 5c\5m, = 5c). (16) 

Note that in the case of general filters, we cannot reduce this 
probabiUty to a half because of the correlation between the 
Fourier modes. Here we define the normalized density con- 
trasts, 



(17) 



(18) 



where the standard deviations Gr and ctq are the root mean 
squares of the variances (^m,)'^'^ {^MiY^^^ which are de- 
fined in Appendix A. Hereafter we refer vidues concerning 
to the subscript r. 

Generally bivariate Gaussian distribution is given by 



1 



27rvT 



:exp 



{vi - ev2f v\ 
~ 2(1 -e2) 2 



di/idt-i, 



(19) 

where the correlation coefficient e is equal to {v\V2) = 
^o(0)/crra-o, which is defined in Appendix A. The conditional 
probability distribution function is obtained by the Bayes' the- 
orem, 



p{V\\v2 = l'2c)^y\ ■■ 



Pi.Vi,V2c) 



1 



:exp 



pilule) 
{v\-ev2cf 
2(1 -e2) 



di^i, 



(20) 



where f2c is the normalized critical density contrast with the 
smoothing scale M2, 5c/<jq, and p{v) is a one-point Gaussian 
normal distribution, 



p{i')dv = ^ exp 



27r 



dv. 



(21) 



In the case of the sharp A;-space filter, the covariance ^o(O) is 
reduced to ct,,, then we obtain e = a^jar. Thus we obtain the 
eq.(15). However in general cases, we must calculate ^o('") for 
specified filters. 

In Figs.l and 2, we show P{M\\M2) = p(i'\\i^2 = i^2c)di^i 
in the cases of the Gaussian filter with the spectral index n = 0, 
and the top-hat filter with n = -2 because e accidentally be- 
comes ao/cr in n = and then P(Mi IM2) = 1 /2. It is evident that 
the Fourier mode correlation increases the value of PiMi IM2) 



above a half when Mi > M*, in contrast to the case of the sharp 
k-space filter, in which the probabiUty is a half for all region 
withM2 > Ml. 



P(Mi\M2) 




Log[M, /M. ] 



6 Log[M,/M] 



Fig. 1. — Probability P{Mi\M2) for the Gaussian filter in the case of the 
spectral index n = 0. The surface of the probability is plotted only against the 
region Mi < M2. 



P(M, 




Log[Mi/M.] 



Fig. 2. — Same as Fig. 1, but for the top-hat filter and for n = -2. 



Solving eq.(14) by substituting P(Mi\M2), we obtain the 
mass functions. In Fig. 3, we show the mass functions for the 
Gaussian and the top-hat filters and the PS mass function in a 
form of a differential multiplicity function, F(M) = Mn(M) / p. 
Note that the function with the sharp fc-space filter is the same 
as the PS mass function. Because P{Mi IM2) increases above 
a half at larger Mi scale while /(> 6c;Mi) is not changed by 
filters, the number density decreases at larger scale. Then, so 
as to cancel out the decrease, the number density increases at 
smaller M2 scale. As a whole, the mass function is seemed to be 
moved to lower mass scale. This result is similar to that given 
by the Monte Carlo calculation by PH90 and BCEK. Thus we 
conclude that the filtering effect can be interpreted as the cor- 
relation effect between different Fourier modes S^. It should be 
noted that in the case of the top-hat filter and n = 0, the correla- 
tion ^o(O) becomes CTq, so the resultant mass function is exactiy 
the same as the PS mass function. 
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(a) n=0 




O 



(b) n = -2 




- — — - Gaussian 
Top— hat 



-6 



-4 



-2 



Log[M/Mj 



Fig. 3. — Differential multiplicity functions F(M) with the filtering effect, 
(a) n = 0. (b) n = -2. The solid, the dashed, and the dash-dotted lines denote 
F(M) for the PS, the Gaussian filter, and the top-hat filter. Note that F(M) for 
the sharp <:-space filter and for the top-hat filter with « = are identical with 
the PS. 



5.2. Peak effect 

Next we derive the mass function under the assumption that 
objects collapse around density maxima. Because ten variables 
are required to describe one peak, we need to calculate eleven- 
variate Gaussian distribution function (one is 6mi and others are 
concerning an object M2). Here we define the first and the sec- 
ond derivatives of the density contrast 5, 

dSix) 



dxi 



(22) 



(23) 



dxidxj 

The peak condition is described by using 77, and Qj as follows: 
T7 = and the eigenvalues of the Hesse matrix <^ are less than 0. 
The covariances concerning 77, and Qj are estimated as follows: 

(24) 



iViVj) = -jSij, 



{x'vi) = 7 = 



(J 



(25) 



(26) 



where x! is a normalized trace of the diagonalized Hesse matrix, 



x' = 



and the variances and correlation are defined in Appendix A. 
Note that (-Ai,-A2,-A3) are the eigenvalues of (J, and other 
three components are transformed to variables given by the Eu- 
ler angles for the diagonalization. 

By using the above quantities, we can calculate the probabil- 
ity P{Mi IM2) as follows, 

P{Mi IM2) = p{vi >Uc\v2 = l^lcV = 0, A;=i,2,3 > 0) 



1-72 



' 27r(l - e2-/i2-72 + 2e^7) 
J^^dx'mi^^di^^cxpi-Q^) 
/o°°jy/(x')exp(-f) 



(28) 



where Qa,Qb and fix) are defined in Appendix A. Substituting 
the above probabiHty to eq.(15), we obtain mass functions. 

It should be noted that in the case of the top-hat filter, the 
second derivative cannot be done because the top-hat filter does 
not smooth out the density field adequately. Then the variance 
for the derivatives of the density contrast diverges to infinity as 
shown in Appendix B. This fact leads that the resultant mass 
function taking into account the peak condition becomes the 
same as that not taking into account it, if we simply interpret 
this divergence as 7, ^ in eqs.(25) and (26). The peak for- 
maUsm should be modified to be able to treat the top-hat filter 
but it is beyond the scope of this paper. 

Apparently the similar thing to the above occurs for the sharp 
^-space filter. If we do not consider the spatial correlation of the 
density field, the peak condition does not add any new proper- 
ties because the behavior of the Markovian random walk owing 
to adding a new Fourier mode conserves. Mathematically this 
is proved by the fact that /x = £7 in the case of the sharp k-spacc 
filter. Thus we obtain the same mass function as the PS even 
taking into account the peak condition. Of course, when we 
consider the spatial correlation, the peak condition will affect 
the resultant mass function because /i does not coincide with 
£7 at r y (see the next subsection). 

In Fig.4, we show PiM\ IM2) in the case of the Gaussian fil- 
ter with the spectral index « = 0. Note that the viewing angle 
is changed. We see that P(Mi IM2) shows a behavior such as a 
function of (M1/M2). Thus the peak condition with the Fourier 
mode correlation increases the value of P(Mi IM2) above a half 
at M2 ^ Ml even at Mi <C M,, in contrast to the results in the 
previous subsection. 



P(Mi\M2) 




-2 -3 
Log[Mi/M.] 



Log[M2/M.] 



(72 



Fig. 4. — Probability P(Mi IM2) with the peak effect for the Gaussian filter 
in the case of n = 0. The viewing angle is changed in order to easily see the 
shape of the probability from Figs. 1 and 2. 

In Fig. 5, we show the resultant mass functions for the Gaus- 
sian filter. For the references, we also show the mass functions 
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of the PS and the top-hat filter derived in the previous subsec- 
tion. In the case of n = 0, the Gaussian filter affects the mass 
function especially at lower mass scale, M < M* and the slope 
at that scale becomes shallower than that of the PS. This change 
reflects the shape of P(Mi jMa) near Mi ^ M2. In the case of 
n = -2, similar to the case of n = 0, the slope at lower scale is 
changed. However the difference from the PS decreases com- 
pared with the n = case. This will lead that the filtering and the 
peak effects become small in the realistic CDM power spectrum 
with n < -2 at galactic scale. Note that PH90 also obtained a 
similar mass function taking into account the peak condition 
with another approximation scheme. 



(a) n=0 




PS 

Gaussian 



o 



(b) n = -2 




PS 

^ Gaussian 

(Top-hat) 

Ln 

I I I I I I I ^ 

-6 -4-2 2 

Log[M/Mj 

Fig. 5. — Differential multiplicity functions F(M) with the peak effect, (a) 
n = 0. (b) n = -2. The solid and the dashed lines denote F(M) for the PS 
and the Gaussian filter. For the reference, the function for the top-hat filter 
only with the filtering effect is shown by the dash-dotted line because the peak 
effect cannot be dealt correctly. As mentioned in Fig.3, F{M) for the sharp 
fe-space filter and for the top-hat filter with n = are identical with the PS. 



5.3. spatial correlation 

Finally, we evaluate the finite-size effect, that is, the effect of 
the spatial correlation on the mass function. As already shown 
in YNG, this effect is similar to changing the critical density 6c 
in the PS to a smaller value than 1 .69, in the case of the sharp 
A:-space filter. Here we generalize the YNG result to general 
filters. 

First we consider only the spatial correlation without the 
peak effect. In this case, the correlation coefficient e must be 
generaUzed to that between two distant points with different 
mass scales, that is, e e(r), and the probability eq.(16) is in- 
terpreted as P(r,Mi|M2). Then by spatial averaging over the 



region of the object M2, we obtain the probability P(Mi IM2), 
J^' P(r,Mi\M2)4Trr^dr 



P(Ml|M2): 



lo'^nr^dr 



(29) 



It should be noted that it is not trivial that the window function 
Wnfiir) is not needed in this averaging. However, since the av- 
eraging operation has been already executed in the calculations 
of and Sm2, we estimate P(M[\M2) by eq.(29). When it is 
included, the effect of the spatial correlation will be enhanced 
because more distant point from the center of the object M2 is 
considered. 




Log[Mi /M, 



Log[M2/M] 



Fig. 6. — ProbabiUty P(.Mi IM2) with the effect of the spatial correlation for 
the sharp ^-space filter in the case of n = 0. 




Log[Mi/M.] 2 ^.g 

Fig. 7. — Same as Fig.6, but for the top-hat filter. 



P{Mi IM2) 



6 Log[M2/M.] 




Log[M,/M.] 



Fig. 8. — Same as Fig.6, but for the Gaussian filter. 

In Figs.6 ^8, we show the spatially averaged probabiUties 
P(Mi IM2) for the sharp fc-space, the top-hat, and the Gaussian 
filters with the spectral index n = 0. In the cases of the two 
former filters, in contrast to the filtering effect, the probability 
becomes less than a half when Mi > M* . The reason why the 
value becomes less than a half is that the correlation between 
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(5mi and becomes weaker at more distant point from the 
center of the object M2. Thus, while the interval of the inte- 
gration of 5mi - 5c is [0, 00] and then the probability is a half 
when the sharp ^-space filter is used and the spatial correla- 
tion is neglected [eq.(15)], we must integrate from larger than 
to cxD under the consideration of the finite volume. Then the 
probability becomes less than a half. On the other hand, in the 
Gaussian filter, the decrease of PiMi IM2) at M\ > emerges 
only at Mi ~ M2. Except for such region Mi ~ M2 > M*, the 
shape of P{M\ [M2) is similar to Fig. 1 in which the spatial corre- 
lation is not taken into account. This suggests that the filtering 
effect is stronger than the effect of the spatial correlation in the 
Gaussian filter 

We show the mass functions for three filters in Fig. 9. While 
the mass functions with the top-hat and the sharp ^-space filters 
horizontally move to larger mass scale compared with the PS 
mass function, the mass function with the Gaussian filter moves 
to smaller mass scale. These properties can be explained from 
the shape of P(Mi IM2). In the sharp A;-space and the top-hat 
filters, P(Mi IM2) is less than a half at Mi > M*. This leads that 
the multipUcity function F(M) becomes larger than the PS at 
that scale because /(> (5e;Mi) in eq.(14) is unchanged. Then 
F(M) decreases at smaller scale Mi < M* so as to cancel the 
increase of F(M) at the larger scale. On the contrary, in the 
Gaussian filter, PiMi IM2) increases above a half in almost the 
region of Mi > M«. This leads that F(M) becomes smaller at 
that scale and larger at smaller scale. 




-6 -4-2 2 

Log[M/Mj 

Fig. 9. — Differential multiplicity functions F(M) with the effect of the spa- 
tial correlation, (a) ra = 0. (b) n = -2. The soUd, the dotted, the dashed and the 
dash-dotted lines denote F(M) for the PS, the sharp i-space filter, the Gaussian 
filter and the top-hat filter. 

Finally we include the peak condition in the above formal- 



ism. In this case, not only e(r) but also /i must be generalized 
to /u(r) because /x(r) describes the correlation between the sec- 
ond derivative of the density maxima with smoothing scale M2 
and the density height with scale Mi distant r from the center of 
the object M2. In order to obtain the probability P(Mi IM2), we 
substitute e(r) and /x(r) into eq.(28) and average over the region 
of the object M2 according to eq.(29). 



P(M, 




Log[Mi/M.] 



Fig. 10. — Probability P(Mi IM2) with the effects of the peak and the spatial 
correlation for the Gaussian filter in the case of n = 0. 



P(Mi\M2) 




-5 .4 
Log[M, /M. ] 

Fig. 1 1 . — Same as Fig. 10, but for the sharp /:-space filter. 



Log[M2/M.] 



In Fig. 10, we show P(Mi\M2) for the Gaussian filter with 
the peak and spatial correlation effects. This figure shows the 
following three characteristics: increase atM2 S> Mi ^M*, in- 
crease toward to M2 ^ Mi , and steep decrease at M2 ~ Mi . It 
can be easily interpreted from Figs.l and 4 that the first char- 
acteristic reflects the filtering effect and the second the peak 
effect. On the other hand, the third is a new property gener- 
ated by the mixture of the filtering, peak and spatial correlation 
effects, because the spatial correlation itself affects P(Mi IM2) 
only at Mi > M^, as shown in Fig.8. Fig. 11 shows P(Mi \M2) 
for the sharp A:-space filter. Similar to the third effect in the 
Gaussian filter case, f (Mi IM2) decreases at M2 ~ Mi . This is 
also a new characteristic. As seen in the peak effect (Fig. 4), the 
change of P(Mi IM2) from a half at M2 ~ Mi affects the slope 
of the low-mass tail of the multiplicity function. 

In Fig. 12, we show the mass functions taking into account 
the effects of filtering, peak, and spatial correlation. For the 
reference, we also show the mass functions for the top-hat fil- 
ter with the effect of the spatial correlation (Fig. 9), because we 
cannot derive those with the peak effect as mentioned in §§5.2. 
At larger scale Mi > M*, the effect of the spatial correlation 
is similar to that without the peak effect shown in Fig. 9. At 
smaller scale Mi < M,, as shown by YNG, the slope becomes 
steeper than that of the PS in the case of the sharp A:- space filter. 
In contrast, in the case of the Gaussian filter, the slope becomes 
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much shallower than the PS. The effect of the spatial correla- 
tion is too weak to cancel out the peak effect which makes such 
a shallow slope, while at larger scale the mass function is close 
to the PS. 

These results show that the slope at smaller scale M\ <C 
is determined by the behavior of P{Mi\M2) from the region 
Ml <i; M2 at which P{M\\M2) is a half toward M\ ^ M2. If 
P{Mi\M2) has a region P{Mi\M2) > 1/2, the slope becomes 
shallow and if P(Mi IM2) < 1/2 near Mi ~ M2, the slope be- 
comes steep. In the case of the top-hat filter, in which the peak 
effect cannot be considered, P(Mi IM2) ~ 1 /2 at M2 ~ Mi < M* , 
then the slope is close to that of the PS. Note that the difference 
of the slope from the PS is prominent in large spectral index, n. 
We will comment on the power spectrum in later section. 
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Fig. 12. — Differential multiplicity functions F(M) with tlie effects of the 
peak and the spatial correlation, (a) n = 0. (b) n = -2. The solid, the dotted 
and the dashed lines denote for the PS, the sharp fc-space filter and the 
Gaussian filter. For the reference the functions for the top-hat filter only with 
the effect of the spatial correlation are shown by the dash-dotted lines. 

5.4. CDM power spectrum 

In this subsection we show the mass function in the case of 
a flat ACDM model as a example of realistic cases. We use a 
power spectrum given by BBKS with fio = 0.3, fi a =0.7 and 
r = 0.21, where JIojOa and F are the cosmological density 
parameter, the cosmological constant and the shape parame- 
ter, respectively. Here we assign the smoothing scale Rq for 
the Gaussian filter to 0.64/? according to BBKS and define the 
mass of objects M = bip-pk'^ according to LC for the sharp k- 
space filter. The power spectrum is normalized so as to be unity 
at 8/i~^Mpc with the top-hat filter. In Fig. 13, we show the mass 
functions for the three filter functions with the effect of the spa- 



tial correlation and the PS mass function. We neglect the peak 
effect in order to compare the mass functions under the same 
condition. 

The behavior for filter changing is similar to the scale-free 
power spectra (see Fig. 9) and the difference among them be- 
comes a little larger. Note that there is an uncertainty in the 
estimation of mass of objects in the cases of the Gaussian and 
the sharp A:-space filters and it cannot be canceled in the case 
of such realistic power spectrum with a characteristic scale. If 
we adopt the same definition of mass as in the previous, the 
difference becomes still larger. 

We also plot the mass function given by Jenkins et al. (2001), 
which is well fitted by their high resolution A^-body simulation. 
At larger mass scale M > lO'^M©, it is in good agreement with 
the mass function with the top-hat filter. This agreement sug- 
gests that the spherical collapse approximation is good in such 
large scale (M > lO'^^Mr;,) and that an overlapping effect be- 
tween similar mass halos becomes important in the intermedi- 
ate scale (M ~ 10''*Mq) as shown by YNG. 




Log[M/h IO'^Mq] 

Fig. 13. — Differential multiplicity functions F{M) with the effect of the 
spatial correlation in a ACDM model. The solid, the dotted, the dashed and 
the dash-dotted lines denote F(M) for the PS, the sharp i-space, the Gaussian 
and the top-hat filters. The crosses indicate the mass function given by Jenkins 
et al. (2001), which was obtained by fitting their iV-body simulation. 

6. DISCUSSION 

In this section we discuss some uncertainties in calculation 
of the mass function. 

6.1. Mass-smoothing scale relation 

So far we use a simple mass-smoothing scale relation, M = 
47r/?^/o/3 and kcR = tt in the case of scale-free power spectra. 
By changing this relation, the mass function may be affected as 
follows. It should be noted that we calculate the mass function 
against the normalized mass M/M* and M* is defined so as to 
cro(M*) = 1 with scale-free power spectra. Thus the mass func- 
tion does not scale in mass by changing the relation, even in the 
case of the Gaussian filter. This effect appears via the spatial 
correlation in the case of the sharp fc-space filter. 

The ratio of the cut-off wave number for the simple rela- 
tion {kc) to that given by the LC method (fec.Lc), kc/kc^LC = 
{2iT^/9)^l^ ^ 1.3. Thus the effect of the spatial correlation is 
weakened because the correlation at distance r is described by 
the combination of kcr, and the probability P{M\ IM2) is ob- 
tained by integrating P(r,Mi|M2) from to /? respect to r. We 
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checked this uncertainty, but the effect does not change our re- 
sults qualitatively, especially at larger mass scale (M > M,). 

If we use a realistic power spectrum, for example CDM mod- 
els, the normalization of the mass spectrum ao(M) becomes not 
arbitrary. In observation the normalization of the power spec- 
trum is given by the amphtude at 8/i"^Mpc with the top-hat fil- 
ter with the bias parameter. This leads that the amplitude at that 
scale does not become unity in cases of other filters. In Fig. 13, 
we adopt the LC method as a usual manner. We have found 
that when we use the simple relation the differences between 
the PS and the mass functions with the Gaussian and the sharp 
A:-space filters becomes large compared with those in the case 
of the LC method. Part of the remaining differences is owing to 
the difference of the amplitude of the power spectrum. 

6.2. Effect of non-spherical collapse 

Throughout this paper, we analyzed the mass function un- 
der the simplification of the spherical collapse, in which the 
critical density contrast for collapse is 5c = 1.69. However, 
in more realistic condition, objects collapse nonspherically. 
Monaco (1995; 1997a, b; 1998) showed that the. first collapse, 
which means one-dimensional or pancake collapse, occurs dur- 
ing more short time-scale. Fitting the PS mass function to the 
resultant mass function, he obtained an effective density con- 
trast 5c — 1-5. This predicts larger number density of halos at 
larger scale compared to the PS, such as the effect of the spatial 
correlation. Apparently this seems to enhance the difference 
of the mass function including the effects of the filtering, peak, 
and spatial correlation. It will be correct when we consider the 
subsequent evolution of baryonic component included in halos 
because shock heating of gas occurs at the first collapse. Even 
so, in analyses of the mass function by using A^-body simula- 
tions which is in good agreement with the theoretical PS mass 
function, it is not trivial that the condition of the identification 
of isolated halos is corresponding to the first collapse. Actually 
such identified halos are rather spherical. This means that the 
effective critical density is larger than 1.5 and that the second 
or third collapse criterion will be better than the first collapse 
one. If so, the non-spherical effect will hardly affect our results, 
unless the mixture of the effects emerges as shown in Fig. 12. 

7. SUMMARY AND CONCLUSIONS 

In this paper, we analyzed the following three effects on the 
mass function: the filtering effect of the density fluctuation 
field, the peak ansatz in which objects collapse around density 
maxima, and the spatial correlation of density contrasts within a 
region of a collapsing halo because objects have non-zero finite 
volume. It has been shown by many authors that the PS mass 
function agrees with A^-body results well. However it had been 
also unclear why the PS mass function works well. Taking into 
account the above three effects and using an integral formula- 
tion proposed by Jedamzik (1995), we tried to solve the missing 
link in the PS formalism. While YNG showed that the effect of 
the spatial correlation alters the original PS mass function, in 
this paper, we showed that the filtering effect almost cancels 
out the effect of the spatial correlation and that the original PS 
is almost recovered by combining these two effects, particularly 
in the case of the top-hat filter As we shown in Fig. 13, the mass 
function with this filter is in good agreement with that given by 
a A^-body simulation (Jenkins et al. 2001). On the other hand, 
as for the peak effect, the resultant mass functions are changed 
dramaticaUy and not in agreement especially at low-mass tails 



in the cases of the Gaussian and the sharp ^-space filters. 

In the sharp fe-space filters, as shown by PH90 and BCEK, a 
density contrast at a point changes as the smoothing scale de- 
creases because a new Fourier mode with an uncorrelated phase 
with other modes gives the density contrast an additional value 
such as Markovian random walk, which originates the so-called 
fudge -factor-of-two. On the other hand, when using general fil- 
ters, we cannot expect that the trajectory of density contrast 
against the smoothing scale behaves like Markovian random 
walk. Then the factor of two is lost. We confirmed the result 
by PH90 and BCEK by using the other formahsm proposed by 
Jedamzik (1995), against the filtering effect which we referred 
in this paper The number density of halos decreases at larger 
mass scale compared to a typical mass scale which gives unity 
variance, and increase at smaller mass scale. Besides we found 
that the slope at smaller scale does not change. Thus apparently 
the shape of the multiplicity function F(M) scales horizontaUy 
against the mass scale. This effect can be interpreted as the 
correlation of the Fourier modes, (5k, in a statistical sense. 

We also investigated the effect of the peak ansatz in which 
objects coUapse around density maxima and which is often con- 
sidered in structure formation theories. In the case of the Gaus- 
sian filter, this makes the low-mass tail of the multiplicity func- 
tion shallower than the PS dramatically. We found that this 
difference can be interpreted by the behavior of the probability 
P(Mi\M2) near Mi ~ M2. P(Mi\M2) increases above a half at 
that scale while it remains a half in the sharp k-space filter. On 
the other hand, in the case of the top-hat filter, we cannot deal 
with this effect correctly because this filter does not smooth 
out the density field adequately so as to be able to differenti- 
ate the density contrasts respect to the spatial coordinate due to 
the sharpness on the edge of the filter. While this trouble may 
be removed to adopt another filter with a rounded edge, it will 
originate complicated mathematical difficulty. 

The effect of the spatial correlation only in the case of 
the sharp A:-space filter was already considered in YNG. They 
showed that, in contrast to the filtering effect shown here, the 
mass function predicts more number of halos at larger mass 
scale, M > M*, and less number at smaller scale, M < M*. 
This is apparently corresponding to scaling the mass scale to 
larger one. This property is caused by decreasing correlation 
between density contrasts with different smoothing scales. By 
considering the other filters, we showed that such properties are 
weakened and the resultant mass functions become close to the 
original PS because the filtering effect strengths the mode cor- 
relation as mentioned above. We also showed, however, that the 
spatial correlation cannot correct the change of the slope of the 
low-mass tail in the case of the Gaussian filter which is caused 
by the peak effect and besides steepens the tail in the case of the 
sharp fe-space filter. This is because the probability f (Mi IM2) 
is mainly affected at Mi ^ M2 almost independent of the ab- 
solute values of Mi and M2 by the peak effect, while mainly 
affected at Mi > M* by the filtering effect and the spatial corre- 
lation. Probably we can confirm which filter is better by using a 
high resolution A'-body simulation, while it may be difficult to 
obtain an adequate dynamic mass range to much smaller scale 
than M* . 

As a conclusion, the discrepancy with the PS shown by YNG, 
in which only the sharp fe-space filter was considered, is almost 
canceled out by the filtering effect. The peak ansatz, which does 
not affect the mass function in the case of the top-hat filter due 
to the mathematical reason, causes a new discrepancy at low 
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mass scale in the cases of the Gaussian and the sharp A;-space 
filters. A new filter similar to the top-hat filter but differentiable 
will be required in order to consider more reaUstic situation of 
collapsing halos. 

It should be noted that the discrepancy is not perfectly can- 
celed by the filtering effect and the spatial correlation, espe- 
cially at larger mass scale, M > M* . This scale is corresponding 
to the galaxy clusters. Several authors have tried to determine 
the cosmological parameters and the normalization factor of the 
power spectrum by comparing the PS mass function to the ob- 
servational cluster abundance (e.g., Bahcall 2000). We should 
pay attention to the fact that the results by such approach may 
depend on the model mass function. 

We discussed an important effect of the non-spherical col- 
lapse, which has been investigated by Monaco (1995; 1997a, b; 
1998). He found that if we adopt the first shell crossing as the 
collapse criterion, the resultant mass function predicts larger 
number density of halos at larger scale M > and smaller 
number density at smaller scale M < compared to the PS 
because the time-scale of collapse decreases. This effect will 
enhance the difference between the PS and the mass function 
given by a full statistics with the effects taken into account in 
this paper and with the non-spherical collapse. Thus when we 
consider the evolution of the baryonic component, we should 
pay attention this effect. However, the agreement of the PS 
with the mass function given by many A^-body simulations will 
be obtained by adopting the second or the third shell crossing as 
the collapse because identified objects as halos are rather spher- 
ical, not sheet-like or filamentary. 

We found that the effects we considered becomes more 
prominent in larger spectral index n. In CDM model, the effec- 



tive spectral index n is less than -2 at galactic scale. Therefore 
in a realistic power spectrum the difference of the low-mass tail 
of the mass function from the PS may be negligible. However, 
at larger scale M > M*, the difference remains. This may leads 
a systematic error when the cosmological parameters are es- 
timated by comparing observations such as cluster abundance 
with the PS. Thus we need more detailed study on the mass 
function. 

Finally we should note that the approximation method used 
here is different from the excursion set approach. In the cal- 
culation of P(Mi IM2), we considered only two different mass 
scales. Ml and M2. However, more strictly speaking, the proba- 
bility must be subject to an additional condition, that is, the ob- 
ject with M2 is an isolated object. In other words, the collapsing 
condition of objects M2 must describe the first upcrossing above 
Sc. This is essentially the same as the cloud-in-cloud problem 
and the overlapping problem which was pointed out by YNG. 
Even so, we believe that the approximation used here is ade- 
quate to resolve the mystery in the PS formalism because the 
fact that the filtering effect correlates different Fourier modes 
leads to enhance the probabiUty that the objects M2 in our for- 
mulation are isolated. This conjecture can be stressed by the 
similarity of the mass functions shown here and the Monte 
Carlo approach to solve the excursion set formahsm with the 
general filters, as mentioned above. In future, this will be 
proved by using a high resolution AT-body simulation numer- 
ical experimentally. 

I wish to thank N. Gouda, T. Yano, and F. Takahara for use- 
ful suggestions. Numerical computation in this work was partly 
carried out at the Yukawa Institute Computer FaciUty. 



APPENDIX 

FORMULATION OF p(mi \m2) FOR GENERAL FILTERS 

In this paper, we only considered Gaussian random fields. It is well-known that the Gaussian field is perfectly determined by 
the two-point correlation function. In the foUowings, we show the form of the probability required to analyze the mass function 
firstly, and show the correlation and the auto-correlation functions, which are introduced as covariances in multi-variate Gaussian 
distribution functions, in the case of the power-law power spectrum, P(k) = (|<5kp) =Ak", where A is the amplitude, which is 
normalized so as to be unity at M = , and the angle brackets stand for the ensemble average. 

The kernel in eq.(14), f (Mi IM2), is generally expressed as a conditional probabiUty. Since we need at least two Gaussian fields 
with smoothing scale Mi and M2 (Mi < M2), or R\ and R2 (R\ < ^^2) in terms of length as mentioned in §2, multi-variate Gaussian 
distribution functions are required. Here we define V'*' = (xi , • • • ,xn) as a vector of A' Gaussian distributed variables and C'^ as a 
condition for A' variables. The distribution function for a variable xn under the condition C'^"^ is written by Bayes' theorem, 

p(x.\C^-')dx,= P^^^^^, (Al) 

where p{Y) is a multi-variate Gaussian distribution function. In this paper, the condition C'^ means a criterion of collapse for larger 
objects, M2. When we do not consider the peak condition, the vector V is two- variable one, = (62, 5i) and the condition C 
constrains one variable 62, C = (^2 = 62c)- Note that the points for S\ and 62 are generally not the same if the spatial correlation is 
taken into account. When we consider the peak condition, eleven-variate Gaussian distribution is required as V'' = {52,'nhCiji5\)^ 
where ; = 1,2,3 and / < j. Note that 77, and Qj are derivatives of 62 respect to the coordinate x,. Conditions for ry, and Qj are 
determined such that the density has a local maximum at the point M2. 

By using a covariance matrix M, the A^-variate Guassian distribution function is generally written as 

P(V^)JVV = (A2) 

V'(27T)'Vdet(M) 

where 

Q = VM-iy'T (A3) 
Mij = {iXi-{xi))ixj-{xj))). (A4) 
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This shows that the muhi-variate Gaussian distribution is perfectly determined by correlations for any two variables. Thus what we 
do is to calculate the covariance matrix M for given variables. 

As shown in YNG, the following covariances have non-zero values: 



6 



(CvCh) = -7\{5ij5ki + 5ik5ji + 5ii6jk), (SiOj) =—r^ij, i'^iCu) =~T^'j' 



15 

where 6ij is the Kronecker's delta. Variances and correlations (7r,<7i, and ^, are defined as 



4tt r 



W(kRi)W(kR2)P(k)- 



sinkr 



kr 



k'^-^dk: 



4ttA 1 
(27r)3 /?f 2'+3 



In+il+iir), 



(27r)3 



AttA 1 

W\kR2)Pik)k"^^dk= ^^—^I„^2M, 



Air f°° AttA 1 

<Tf = TT— 3 / W\kRi)P{k)k^dk = ^In^l. 

(.^""J Jo 



(A5) 

(A6) 
(A7) 
(A8) 



/o (27r)3 /?f 3 

Clearly ^/(r) means the 2/-th moment of the two-point correlation function, is the 2Z-th moment of the variance with the smoothing 
scale /?2, and al is the variance with the scale Ri. The subscript r stands for the value with the scale Ri because generally it is 
separated by r from the center of the object with the scale R2 in this paper. It should be noted that ^/C) — * o-^ for r — > and R\ — > R2. 
The integrals /,„(r) and /„, are calculated for specific window function WikR) in the next section. 

In the followings, the peak condition is taken into account. In the case without the condition, the probabihties are obtained by 
neglecting quantities concerning 77, and 0;- Here we diagonalize the Hesse matrix Qj by introducing the eigenvalues (-Ai,-A2,-A3) 
with the order Ai > A2 > A3, and the Euler angles a', /?' and 7'. The volume element concerning the Hesse matrix becomes 



W d0j=\iXi-X2)iX2-X3)ih-M)\dXidX2dX3smf3'dP'da'dy. 



ij= 1,2,3 



Here we normahze the variables. 



Si 62 , A1 + A2 + A3 , A1-A3 , A1-2A2+A3 

Vl = —,U2 = —,X = ,3' = ,Z = . 

Ur (To (72 ZCr2 ^CT2 



(A9) 



(AlO) 



The peak condition constrains the derivatives of (5 to be jy = and each diagonal component —A,- is less than 0, which means —y' < 

z' < y' and y' < 0. The covariances are transformed as 



1 



(1.2) = 1, = 1, (XV2) = 7, (^^11^2) = eir), (xVi) = nir), {x'^) = 1, (y'^) = —,{z'^) = - 



15' 



where e(r), /x(r) and 7 are written as 



e(r) = {U1U2} = 
^l{r) = {x'vl) 



7 = {x'v2) = 



6W 
6(0 

(Tr<T2 ' 
O'ocr2 ' 



and /u(r) ^ 7 for r ^ and Ri -^R2. In order to obtain the exponents of probabilities, we define the following variables: 



Qa = 
Qb = 



1-7^ 



1 - - /Lt^ - 7^ + 2e/X7 

ix'-7l^2? 



2^2 + -, 2^ 

1-7^ 1-7^ 



1-72 



(All) 

(A12) 
(A13) 
(A14) 

(A15) 
(A16) 
(A17) 



e, = zv| + 15/^+5z'^ + 3^, 

where we omit the argument r of e(r) and nir) for simplicity. After integrating the volume element over the Euler angles, we obtain 
p{i'i,i'2,r),x' ,y' ,z')di'idi'2d^r]dx'dy'dz' 

Qa + Qb + Q. 



3(15)V2|y(v'2_-'2^| 



8 V27rV2(T3 i/l-e2-y^2_^2 + 2e;Li7 



p{v2,ri,x' ,y' ,z!)dv2d^ridx'dy'dz! 



3(i5)^/2|y(y2-z^2)| 

87r3 (73^1-72 



exp 



Qb + Qc 



exp 



dv2d rjdx'dy'dz' ■ 



dvidv2d^ ridx'dy'dz , 



(A18) 



(A19) 
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The latter equation is the same as eq.(A7) in BBKS. From the Bayes' theorem, the required probability is written under the above 
condition against z/2,T7,x',y and z' as 

Pir,Mi IM2) = J , / -'-^ on ' ^^20) 

y 27r(l-e^-/Lt^-7^ + 2e/i7) dx' f(x')e-Qi>/^ 



where 



fix') = / dy' dz'+ dy' / dz! |/(/2-z'2)|exp • 

\J0 J-y' Jx'/A Jiy'-x' J \ 

1 5x' exp + VTO^ I erf K j + erf L' ^ 



2250 

with the definition of the error function, 



(A21) 



erf(x) = ^ / e-'^dt. (A22) 



Note that the form of f{x) is different from eq.(A15) in BBKS because they considered not the probability but the number density of 
peaks. Finally, if we take into account the spatial correlation, that is, the different points of 61 and 62, we need to spatially average 
the above probabiUty over the region R2 in order to obtain P{M\ IM2), 

P{M,\M2)= ^' '' I' . (A23) 

CORRELATION COEFFICIENTS FOR SPECIFIC FILTERS 

In the foUowings, we give the specific forms of the correlation coefficients e(r),/x(r) and 7 for the three filters, the top-hat, the 
Gaussian, and the sharp ^-space filters, and for the spectral index n = and -2, which are important values when considering 
formation process of galaxies and galaxy clusters in a CDM universe. 

Here, for convenience, we define the following dimensionless variables, 

y = kR2, z = Ri/R2, s = r/R2. (Bl) 
By using these variables, the correlation coefficients are written as follows: 

7= , , e(r)=z2— , ii(r) = z^ ^==. (B2) 

V 'n+2'n+6 V in+2ln+6 

Thus, the problem is reduced to calculating the integrals /m(0 and /„. 

Top-hat filter 

In the case of the top-hat filter, the behavior of the correlations changes at the boundary of an object with the smoothing scale R2. 
So we define the regions 1, 11, and 111 as 

1. Internal: r <R2-R\, 

11. Crossing: R2-R\ < r<R2+Ru 

III. External: R2+R\ < r, 

and the following new variables as 

a=l-z-s, b=l-z+s, c=l+z-s, d=l+z+s. (B3) 
Actually the region III is not considered in our analyses, but we express the expUcit forms of the correlations for the region III in the 
following. 

The generic forms of the integrals can be expressed as 



f°° 9 sinys 
Imir)= / -T^{sinyz-yzcosyz)(siny-ycosy)y"' dy, (B4) 

Jo y ^ 

Im = / — (sin.y->xos3;)V"t^3'- (B5) 
Jo >' 

The solutions for m = 0,2,4 and 6 depend on the region 1, 11 and 111 as follows: 

k(r) = I - 6 1 fl |( 1 - z)a^ + sign(fl)(fl*' - 3Qza*) 

+ 6\b\(l- z)b'^ + sign(b)(-b'' + 30zb^) 
+ 6\c\(l +z)c'*+sign(c)(-c^-30zc'*) 

- 6\d\il+z)d'^+signid)(d^ + 30zd% (B6) 



NAGASfflMA 



/2W = 4|fl|(l-z)a2+sign(fl)(V+12za2) 

-4\b\(l-z)b^ + sign(b)(b*-l 2zb^) 
- 4 1 c |( 1 + z)c2 + sign(c)(c'* + 1 2zc^) 
+ 4 1 (/ 1 ( 1 + z)(f + sign(c/)(-c/'* - 1 2z(f)] , 



or more explicitly, 



97r 

^4('-) = - 2|fl|(l -z) + sign(fl)(fl2-2z) 



+ 2\b\{\-z) + s,ign{b){-b^ + 2z) 
+ 2 1 c I ( 1 + z) + sign(c)(-c^ - 2z) 

-2|^;|(i+z)+sign(^;)(^;2+2z)], 



97r 
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/eCO = [-sign(a) + sign(^) + sign(c) - sign(<i)] , 



(-3z'+15-5/) 

[ 1 0(z*' - 9z'* + 1 6z^ - 9z^ + 1 ) + 48 s(-z^ + 5z^ + 5z^ - 1 ) 



Io{r) = { 640sz^ 

+90s\z^ - 2z^ + 1) + 80s\-z^ - 1) + 30/(z'* + 1) - 2**'] 

TT 



and 



^'^"^^ i ^[-3il-z^f + Ssil+z')-6s\l + ^)+s'] 



hir) = 



9tt 



(z^-s^ + l) 



l2 = 



37r 
T' 



I 

II 

III, 



I 
II 

in, 



II 
III, 



l4= lim — [-2 + 2)' +2cos2)'+ysin2y] ^ oo, 

y— >oo 4y 

4= lim ;^[12);+4);^ + 18jcos2);-15sin2);+6);^sin2);] ^ oo. 

The correlation coefficients are obtained by using the above integrals for the spectral index n = 0, 

7 = 0, 



/X 3/2 ^2(r) 



^3/2 
7-3/2 



lo 



I6s 



(-3(l-zV+8i(l+z^)-6/(l+z2)+/) 



and for n = -2, 
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e(r) = zV2 



37r/5 



^(-3^2+15-5/) 

-5/2 

:r— [ 1 0(z^ - 9z^ + 1 6z^ -9z^+l) + 4Ss(-z^ + 5z^ + 5z^-l) 
384s 



+90i^(z'' - 22^ + 1) - SOs\z^ + 1) + 30/(z^ + 1) - 2s^] 



6s 



n 
III, 



l/2_^W_ 



/V4 



0. 



(B22) 



(B23) 



These show that we cannot deal with the configuration of density peaks, at least statistically, in the case of the top-hat filter because 
of the divergence of I4 and I(, to infinity. This may be originated by the sharpness of the window function, which do not adequately 
smooth out density fields so as to be differentiable. 

Here we show that the peak condition is canceled out in the calculation of the probability P{r,M\ IM2). Now fi{r) = 7 = 0, then the 
eq.(A20) becomes 

1 \r. dx f{x ) dv\e 2 

JO 'Juu ^g24) 



and the exponents become 



P{r,Mi\M2) = 



Qa 



{vi-ev2cf 



l-e2 



Because the exponent Qa does not depend on x!, we can cancel the integration respect to x! as follows. 



P(r,Mi|M2): 



V27r(l-e2) 

This is the same form as the probability not taking into account the peak effect. 



dvie 



(B25) 



(B26) 



Gaussian filter 

The generic forms of the integrals are expressed by using the confluent hypergeometric function and the gamma function as 
follows: 



Specifically, 



and 



ys 



[2(1 



l+m 3 



2 '2' 2(1 



f 

Jo 



fdy. 



loir) = -erf 
2s 



h{r) = 
h{r) = 

hi.r) = 



V2(l+z2) 
1 



r(^) 



2(l + z2)3/2 

TT 3 
2(l + z2)5/2 

TT 15 
2 (l + z2)7/2 



n4 



3(1 

2s^ 

3(1 +z2)"^ 15(1 +z2)2 



e 2(1+^'), 



3V7r 150r 



The correlation coefficients for n = are 



7: 



3 



3/2 



(B27) 
(B28) 

(B29) 

(B30) 
(B31) 

(B32) 

(B33) 
(B34) 
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and for n = -2, 



7 = 
Kr) = 



1 

7!' 



erf 



1/2 



3(l+z2)3/2 



(B37) 

(B38) 

(B39) 



Sharp k-space filter 



Here we show the integrals by using the generic cut-off wavenumber kc which is introduced in the filter, yc = kcR. In this paper we 
showed the results for = tt. If the definition by the LC method is used, jc = (97r/2)^/^. The generic forms of the integrals are 



ys 



1 +m 



l+m 3 3+m (syc)^ 
2 '2' 2 ' 4~ 



-'m Jc I 



where \F2{a;b,c\z) is the generalized hypergeometric function. The specific forms of Imif) are 



1 sin/ , 
- / —dt, 
s Jo t 



1 



h{r) = -7 1 sin .?>•,. - syc cos syj , 



1 



h{r) = -r[syc{(>-s y ^) cos, sy c - 2>{2 - s y^)s,msyc\, 
k{r) = \[-sycil20-20sy,+syjcossyc + 5i24-l2s'yl+sy,)smsy,]. 



The correlation coefficients forn = are 

e(r): 



5 ' 
3z3/2 

-r-r [smsyc - syc cossyc] , 

— [syci6 - s^) cos syc - 3 (2 - />'^) sin syc] , 



and for n = -2, 



7: 



V5 

3 ' 

7I/2 



sm/ , 
df, 

syc Jo t 

, , [smi);^ - yyc cos yyc] • 



(B40) 
(B41) 

(B42) 

(B43) 
(B44) 
(B45) 

(B46) 
(B47) 

(B48) 

(B49) 
(B50) 
(B51) 
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